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Abstract: We develop the helicity formalism for spin-2 particles and apply it to the 
case of gravity in flat extra dimensions. We then implement the large extra dimensions 
scenario of Arkani-Hamed, Dimopoulos and Dvali in the program AMEGIC++, allowing for 
an easy calculation of arbitrary processes involving the emission or exchange of gravitons. 
We complete the set of Feynman rules derived by Han, Lykken and Zhang, and perform 
several consistency checks of our implementation. 
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1. Introduction 

We have learned from string theory that extra spatial dimensions are the price for unifying 
the Standard Model with gravity. Yet we are practically ignorant about how and why these 
extra dimensions are hidden from our world. Will they appear at the electroweak scale? 
Are they bosonic or fermionic? We hope that the answers to these puzzling questions are 
awaiting us in the next round of collider experiments at the Tevatron, LHC, and, possibly, 
at the next linear collider. 

If spatial extra dimensions exist and are accessible to at least some of the Standard 
Model fields, they have to be very small, < C(TeV _1 ). However, it is conceivable that all 
Standard Model particles are confined to a brane, and only gravity (which is rather poorly 
tested at short distances) is allowed to permeate the bulk. This is the classic scenario of 
large extra dimensions [1,2], from now on referred to as ADD. In ADD the experimental 
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limits on the size of the extra dimensions are much weaker. For recent compilations of 
various experimental probes of large extra dimensions, see e.g. [3-5]. 

Such scenarios offer an alternative solution to the gauge hierarchy problem and new 
avenues for unifying the Standard Model interactions with gravity. They also predict a 
variety of novel signals, which are already testable at high-energy colliders, see [6] and 
references therein. Generally speaking, there are two classes of signatures: real graviton 
production [7-12] and virtual graviton exchange [8,9,13]. Gravitons have enormous life- 
time and decay off the brane, hence real graviton production results in generic missing 
energy signatures. Virtual graviton exchange may generate numerous higher dimension 
operators, contributing to the production of fermion [8,9,13,14], gauge boson [15-17] or 
Higgs boson [18, 19] pairs, as well as more complicated final states with more than two 
particles [20,21]. 

The phenomenology of ADD extra dimensions has been studied rather extensively, 
and numerous experimental bounds from LEP and the Tevatron on the scale of extra 
dimensions M s have been set, yet none of the existing general purpose Monte Carlo event 
generators has fully and consistently implemented the ADD model. For example, PYTHIA 
[22] now allows only for the production of narrow graviton resonances as in Randall- 
Sundrum-type models [23,24]. Virtual graviton exchange is not implemented. For Run 
I analyses [25-27] of large extra dimensions, the Tevatron collaborations have been using 
an unofficial version of PYTHIA, where ADD graviton production has been implemented as 
an external process [28]. The situation with HERWIG [29] is very similar - narrow graviton 
resonances are included [30,31], but not much else. ISA JET [32] has recently added the 
process of real ADD graviton production in association with a jet or a photon [33] and 
the implementation is analogous to [28]. PANDORA [34] has implemented virtual graviton 
exchange for e + e _ — > ff,^^,ZZ,W + W~ with the contact interaction of [13] as well as 
e + e _ — ► as in [7]. Finally, programs such as Madgraph [35] or COMPHEP [36], although 
highly automated, do not contain spin-2 particles and thus cannot be applied for graviton 
studies in their present form. Given all this, the pressing need for an automated event 
generator, capable of studying ADD phenomenology in its full generality, becomes quite 
evident. 

The ADD Feynman rules have appeared in [8,9] 1 , so it may seem that implementing 
the ADD model in one of the above-mentioned event generators is rather straightforward. 
However, apart from a few simple cases of two-body final states, the calculations involve 
a large number of diagrams, whose evaluation through the standard textbook methods 
of squaring the amplitudes, employing the completeness relations and taking the traces, 
becomes prohibitively difficult. This problem is considerably alleviated when the helicity 
method is used, which avoids lengthy matrix multiplications stemming from the notorious 
chains of Dirac matrices. First ideas pointing into this direction can be found in [37,38]; 
later on, different approaches were formulated in [39-49] and further refined in [50-52]. 
Thanks to their easy implementation and automatization (see for instance [53-55]) helicity 
methods have become a cornerstone for modern programs dealing with the evaluation of 



1 See, however, section 3.2. 
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cross sections for multi-particle final states. Examples for such programs are the multi- 
purpose matrix element generators GRACE [56], HELAC [57], ' MEGA [58], AMEGIC++ [59], 
ALPGEN [60], and Madgraph [35]; furthermore there is a large number of more specialized 
codes such as those described in [61-63]. 

The main purpose of this paper is to develop the helicity formalism for spin-2 particles 
and to implement the ADD scenario in AMEGIC++ (A Matrix Element Generator In C++). 
Previous studies of ADD extra dimensions using helicity amplitudes can be found in [12,34]. 

The outline of the paper is as follows. In Section |2| we review the basics of the helicity 



formalism and discuss some ideas (2A and ^J5|) that allow a general extension to the case 
of spin-2 particles (gravitons). In Section |3] we specify to the case of ADD and highlight 
some features of the specific implementation in AMEGIC++. For the main part we follow 
the conventions of [9]. We describe both external and virtual graviton processes and list 
a couple of Feynman rules that had not been published before. Section || compares this 
implementation to some existing results in the literature, while Section |B| contains our 
conclusions. Several novel studies of ADD phenomenology at current and future colliders 
will be presented in a separate publication [64]. 



2. The Helicity Method 

The basic idea behind the helicity method is to replace Dirac matrices by suitable spinor 
products, and to use the fact that expressions like 

n(Ai, pi)u(X 2 , V2) 

are simple scalar functions of the helicities A and the momenta p of the spinors. 
2.1 Spinor basis 

The spinor basis used in AMEGIC++ follows closely the definitions of [51,52], where spinors 
are defined in such a way that they satisfy the Dirac equations 

$u(p) = +m ■ u(p) , $v(p) = —m • v(p) , (2-1) 

with 

j2 2 2 
p = p = m , 

where m 2 is not constrained to be positive. For space-like p, i.e. p 2 < 0, m is purely 
imaginary, and the resulting phase degeneracy can be lifted by considering eigenstates of 
75^. With the polarization vector s satisfying 

s-p = 0, s 2 = -l, (2.2) 

this translates into fixing the degeneracy by the conditions 

7sMp> ±) = ±U (.P, ±) > 75^(p, ±) = ±v(p, ±) . (2.3) 
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Spinors can now be constructed explicitly along the lines of [47]. Denning two orthogonal 
auxiliary vectors ko i with 



k 2 = , k\ = -1 , k ■ h = , 
massless spinors w can be defined through 

w(k , X)w(k Q , A) = 1 + 2 ^ 75 $o , 
and the relative phase between spinors with different helicity is given by 

w(ko, A) = X^iw(ko, —A) . 



(2.4) 



(2.5) 



(2.6) 



Spinors for an arbitrary four-momentum p with p 2 = m 2 can be obtained from the massless 
ones with the help of 



u(p, A) = f + w(k , -A) , v(p, A) = -= m w(k , -A) 



V 2 P ■ k o 



V 2 P ■ k 



(2.7) 



both satisfying Eqs. Q2.1| ) and ( ^3| ) 2 . In this approach, conjugate spinors are defined to 
satisfy relations similar to those of Eq. fl2.3|), 



u(p, ±h 5 ^ = ±u(p, ±) , v(p, ±)7 5 ^ = ±«(p, ±) , 
and their relation to the massless spinors is given by 

u(p,X) = w(k ,-\)-^=^= , v(p, A) = w(k , -A)- ^ 



This guarantees the orthogonality of the spinors and after normalizing 

u(p, X)u(p, A) u(p, A)u(p, A) 

2m 2m 

the following completeness relation holds true: 

1 = ^ U (P> X )u( P> A ) ~ V (P> X)v{p, A) 

A 

2.2 Fermion lines 



2m 



(2.9) 



(2.10) 



(2.11) 



In particular, the numerators of fermion propagators experience the following manipulation 



(j, + m ) = ^2 



1 + 



m 



u(p, X)u(p, A) + 1 



m 



v(p, A) v(p, A) 



, (2.12) 



which is valid also for space-like p, i.e. for p 2 < 0. For a more detailed derivation, see [52]. 
This replacement allows us to "cut" diagrams into pieces with no internal fermion lines. 



2 In [65] spinors for massive fermions have been proposed which are directly related to physical polariza- 
tion, i.e spin, states. 
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2.3 External gauge bosons 

For spin-1 bosons there are two different strategies to deal with the polarization vectors: 

1. either replace them by spinor products, or 

2. define explicit polarization 4- vectors. 
Both options are available in AMEGIC++. 

2.3.1 Replacing polarization vectors by spinor products 

In this section, the replacement of polarization vectors of external massless or massive 
spin-1 bosons by spinor products is briefly summarized. 

• Massless spin-1 bosons on their mass-shell have two physical degrees of freedom; 
hence their polarization vectors satisfy 

e^ P ,X)p" = 0, e fl (p,X)e^p,X) =0, 

e>, -A) = e* (p, A) , e„(p, A) d>(p, -A) = -1 , [ ' ' 

where p and A are the four-momentum and the polarization label of the boson, re- 
spectively. In the axial gauge the polarization sum reads 

£ e>, A)e>, A) = + Ma+P^t , (2.14) 

A=± 1 H 

where the auxiliary vector q iff p and r/^ denotes the flat metric tensor. Then, it can 
be shown that the replacement of the polarization vectors through 

fy(P, A) = 1. u(q, A)7 M n(p, A) (2.15) 

y/Ap-q 

satisfies all relations of Eqs. ( |2 . 1 3j ) and ( |2.14| ) and is thus a good choice. 

• For massive vector bosons one chooses a similar strategy to construct a spinor ex- 
pression replacing the polarization vectors. However, massive spin-1 bosons on their 
mass-shell have an additional degree of freedom, the two transverse polarizations are 
supplemented with a longitudinal one with A = 0. This state obeys 

e M (p,0)p" = 0, ^(p,0K(p,A) = 0, e M (p,0)e'*(p J 0) = -l ) (2.16) 
and the corresponding spin sum becomes 



£ e,(p, X)et( P , A) = -n^ + ^ , (2.17) 

A=±,0 



where p and m are the four-momentum and the mass of the boson. In general, there is 
no simple replacement for the polarization vectors as for the massless case. However, 
for unpolarized cross sections, the only relevant requirement is that Eq. (2.17) be 
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reproduced. This can be simply achieved by introducing vectors and kinematics 
related to a pseudo-decay of the boson into two massless spinors with momenta n 2 
satisfying 



r 



1 = 0, r» + r%=p». (2.18) 



Obviously the quantity 



= u(r 2 , -)7 Ai -u(r i , -) (2.19) 



projects out any scalar off-shell polarization. This quantity can be used to replace 
the polarization vector, because the integral over the solid angle f2 of n available for 
the decay leads, up to a constant normalization factor, to the spin sum of Eq. ( 2.17| ), 

^J^a'^-^ + Mf. (2.20) 

Therefore, replacing 

e M ^a M , e » e t^ g— -2 / dVLa^al (2.21) 

A=±,0 J 

leads to the desired result. A first discussion of these ideas can be found in [45]. 
2.3.2 Polarization basis 

It often turns out to be more convenient to keep the polarizations explicitly. A possible 
choice for an on-shell vector boson with momentum 

Vn = (Po> |p| sin 9 cos </?, \p\ sin 8 sin </?, \p\ cos 9) 



is 



e fi (p) = "75 cos ^ cos L P zs r^ srn ^) cos $ s i n ' cos 99, — sin 0) 

e »=^f(b1^o|V (2 - 22) 



For the special choice of q = (po, —p) this satisfies Eq. (2.13) and Eq. fl2.14| ) for a massless 
boson; for massive bosons, also Eq. ( |2.16| ) and Eq. ( p. 17 ) are fulfilled. 

Note that from now on we will label the polarization mode with an upper index. To 
avoid confusion, Lorentz indices are always denoted by /u, z/, p or a, while A and n are used 
to indicate a vector polarization mode. 

2.4 Gauge boson propagators 

Similar to the external vector boson treatment, there are two strategies to deal with gauge 
boson propagators 3 : 



3 This refers to the soon-to-be released version 2.0 of AMEGIC++. 
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Structures in Feynman amplitudes with internal boson lines are calculated as one 
function (composite building block), depending just on the spinors at its edges. (Ex- 



ternal vector bosons are replaced as in 2.3.1 



• The numerator of vector boson propagators can be replaced by products of polariza- 
tion vectors similar to the completeness relations in Eq. ( 2.14j) and Eq. ( [2.17 ). This 
allows us to "cut" internal boson lines. 

The replacement of the numerator for a massive spin-1 boson in unitary gauge reads 



X=tl,t2,l,S 



where p^ = (po, \p\ sin 6 cos (p, \p\ sin 6 sin ip, |p| cos 6) can now be off-shell. The correspond- 
ing polarization vectors are 

ej^(p) = (0, cos 9 cos (p, cos 9 simp, — sin 9) , 
e ^{p) = (0, -sin <^9, cos 99,0) , 




"» = \\- p T^P^ ( 2 - 24 ) 

with the physical boson mass m. To simplify the implementation we choose a basis with 
linear polarization vectors e^ 1 and e^ 2 for the transverse modes. Together with the longi- 
tudinal mode these vectors yield Eq. ( [2.17 ) for particles on their mass-shell. They are 



normalized and orthogonal. The scalar polarization reflects the "off-shellness" of the 
propagating particle. The orthogonality to all other polarization vectors is obvious, since 
it is proportional to p^. Note that in Eq. ( p. 23 ) there is no complex conjugation in the 



polarization product. This avoids an extra sign discussion for space-like p, where and 



may become imaginary. The above replacement is also valid for vector bosons with a 
width, if m 2 is substituted by m 2 + imT. 

A simple inspection of Eq. ( [2.23 ) shows that a similar relation for massless vector 



bosons can be obtained by taking the limit m — > 00. Then, only is affected and turns 
out to be 



tf(p) = y-^JV ( 2 - 25 ) 

By contracting the corresponding Lorentz index in each of the neighbouring vertices with 
a polarization vector and summing over the degrees of freedom, the Lorentz structure of 
the propagator is absorbed. 
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2.5 Spin-2 bosons 

In this section we will show how the above ideas can be extended to deal with spin-2 
particles. We follow the conventions and Feynman rules derived in [9]. With some minimal 
changes, this can be applied to other models involving quantum gravity as well [23,24]. 
The propagator for a KK-graviton with mass m is 

iA W,M(^ ^ . (2-26) 



where 



o , x 1 / P»Pp\ ( PvPu\ . 1 / PnP*\ ( PvPp 



Employing Eq. ( 2.23| ) for each term in parentheses, we obtain 



tfpu,pa ~ 2^ I 2 M P f cr 2 I* a v P ~ ~q (i v e p € a J ■ ( 2 - 28 ) 

X,K=tl,t2,l,S 

Since the graviton couples to the energy-momentum tensor, all resulting vertices are sym- 
metric in /i *-* za Using this symmetry the above equation can be simplified to 

Bpu,pa = ^ ' i e p e p e u e a ~ ~^ e p e v e p e <j J ■ (2.29) 

X,K=tl ,t2,i,S 

The polarization vectors e A have been defined in the previous section. To absorb the 
Lorentz structure of a spin-2 propagator into the vertices, their Lorentz indices now have 
to be contracted by two polarization vectors and the double sum in Eq. ( gjp has to be 
performed. 

External spin-2 particles can be described by polarization tensors. A possible choice 
is given in [9] 4 : 

k = (++ J_ / + , o +\ J_f,+,- , e - e +_2e°e°) —(f-f +e°e-) e~e~\ 

^pV yp-V, V P V < C p C U J 1 /g \ C p C V • C„ ^*ZyX v ) , ^C^ C„ T c^Cj, . 

(2.30) 

Again applying the [i <-> f symmetry in the vertices, this simplifies to 

4, = le+e+, V2e+el , ^| (e+e; - e°e°) , V2e^° u , ^ e ; | . (2.31) 
Here, e^'° are the usual polarization vectors for massive vector bosons, which were defined 



in Section |2.3.2| . All polarization tensors are normalized and orthogonal, and they fulfill 
the completeness relation 

€ pu e pa = B^ U)Pa . (2.32) 



fc=l...E 



4 Notice that our normalization convention differs from the last version of [9] by a factor of y/2. We chose 
to do this in order for the polarization sums to be normalized to 1, as in Eq. (A. 3) in the published version 
of [9]. 
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x 1 x 2 




X\X 2 


Y(pi, Ai;P2,A 2 ;cL,Ci J ) 


++ 


CRHlT]2 + C L fl2m 


+ - 


clS(+;pi,p 2 ) 



Table 1: Y- functions for different helicity combinations. Missing combinations can be obtained 
through the replacements + <->■ — and L <-> i?. 



2.6 Construction of building blocks 

Utilizing the manipulations discussed above, all possible Feynman amplitudes can be built 
up from a small set of basic building blocks. Denoting all spinors by u, denning chirality 
projectors Pl,r through 

1=F75 



PLJt 



and defining 



mi 



Vi 



rji = sj 2pi ■ k , m 

where +/— refers to particle/antiparticle, the basic building blocks are as follows: 
1. Basic functions are direct products of massless spinors of the form 

S(X;pi,p 2 ) = u(pi,A)w(p 2 ,-A) . 
Expressed through the auxiliary vectors ko,i, they read 

±(Pi ■ ko){P2 ■ h) =F (pi • fei)(p2 • feo) - it^p°P\P2K k \ 



(2.33) 
(2.34) 



(2.35) 



^(rh;^!,^) 
5(±;pi,p 2 ) 



r?ir/ 2 



-5(±;p 2 ,Pi) 



(2.36) 



2. The simplest structure is the direct product of two (possibly massive) spinors of the 
form 



Y(px, XDP2, \ 2 ; c R , c L ) = n(pi, Ai) [c R P R + c L P L ] u(p 2 , A 2 ) 



(2.37) 



which emerges for instance through couplings of scalar bosons to a fermion current. 
The results for this function are summarized in Table |]. 

3. The next basic structure is the product of two spinors with a Dirac matrix in between. 
This emerges for instance when external gauge bosons couple to fermion lines with 
the polarization vector of the gauge bosons expressed explicitly or from terms ~ p^p u 
of massive gauge boson propagators. This building block is called X-function and 
reads 



X(jn, Ai;p25P3, A 3 ;c l ,cr) = u(p\, Ai )^ 2 [cl -Pl + crP r ]u(p 3 , A 3 



(2.38) 



Only the restriction kop 2 7^ has to be applied by a suitable choice for k$. The 
results for this function are summarized in Table 0. 
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A1A3 




++ 
+- 


ViV3V2 c L + f4vimcR + c R S(+;pi,p 2 )S{-;p2,P3) 
CLliiV2S(+;p 2 ,P3) + c R n 3 rj 2 S(+;p 1 ,p 2 ) 



Table 2: X-functions for different helicity combinations. Missing combinations can be obtained 
through the replacements + <-> — and L <-> R. 



A1A2A3A4 


Z (pi , Ai ; p 2 , A 2 ; Ps , A 3 ; p 4 , A 4 ; c^ 2 , c} 2 ; c| 4 , c| 4 ) 


+ + ++ 


2 [5 , (+;p 3 ,Pi)5'(-;p2,P4)c) ? 2 c| 4 + m^rj^c^c^ + ^/M^i^e^cf 4 ] 


+ + +- 


2t? 2 4 2 [5 , (+;pi,p 4 )At3c| 4 + s '(+;Pi,P3)At44 4 ] 


+ + -+ 


2r?ic)j 2 [-S(-;p3,P2)M4^ + 5 , (-;p 4 ,P2)/"3c| 4 ] 


+ + — 


2 [S{+;p 4 ,pi)S{-;p2,P3)c x R c^ + m^mm^L + ^3^ViV2C R c| 4 ] 


+ - ++ 


2r ?4 c| 4 [5(+;pi,p 3 )^24 2 + 5(+;p 2 ,P3)/^ici 2 ] 


+ - +- 





+ - -+ 


-2 [/j, 1 fj,ir)2r]3C^(^ : + ^ 2 /W774cJf c| 4 - /i 1/X3 772 774c} 2 c| 4 - /^/t^i^c^cf 4 ] 


+ 


2 %c| 4 [5'(+;P4,P2)Mici 2 + 5"(+;pi,p 4 )/j 2 c)j 2 ] 



Table 3: Z-functions for different helicity combinations. Missing combinations can be obtained 
through the replacements + <-> — and L <-> i?. 



4. The last structure to be discussed is 



Z(pi, \i;p 2 , A 2 ;P3, A 3 ;p4, A 4 ; c} 2 , c|j 2 ; c| 4 , c R ) = 
u( Pl ,\ 1 )^[c?P L + 4 2 Pr]«(p2, A 2 )n(p 3 , A 3 ) 7 , t [c| 4 J PL + cfP R ]u{jp Al A 4 ) ,(2.39) 

which stems, for example, from terms ~ if"' in gauge boson propagators between 
fermion lines. Results for this function can be found in Table H. 



3. Implementation of the ADD Model 

In our implementation of the ADD model we have applied the formalism developed in [9]. 
We assume n flat extra dimensions of common size R, compactified on an n-dimensional 
torus. While the size R is unambiguous, it has become customary in the literature to quote 
bounds on the size of the extra dimensions in terms of an energy scale M, for which there 
are several popular conventions [8,9, 13]. We choose to work with the scale Ms as defined 
in the revised version of [9]: 

(4vr)"/ 2 r(n/2) 

R ~ 2M^G N ' (3 ' 1} 

where Gn = 1/(8ttM p ) is the 4-dimensional Newton's constant, with Mp = 2.4 x 10 18 
GeV the reduced Planck scale. The other commonly used scale, Mp, introduced in [8] is 
related to Ms as 

M s = 2v^F [I>/2)] 1/(n+2) M D . (3.2) 
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Virtual graviton exchange processes are sometimes parametrized in terms of higher dimen- 
sion operators whose coefficients are suppressed by a scale [8] or Mi [13]. The relation 
to Ms can be obtained, comparing ( |3.9| ) with ( 3. 1C| ) and ( |3.11| ). 



3.1 Summation of the KK states 

While the coupling of matter particles to a single graviton is suppressed by the inverse of 
the 4-dimensional Planck mass Mp, the amplitude of any given process involving gravitons 
is enhanced by the multiplicity of graviton states. After summing over the tower of the 
graviton mass states, the effective suppression is only of order M$. 



3.1.1 Virtual KK particles 

The amplitude for any process involving a virtual graviton exchange contains the sum 



rat + ie 



(3.3) 



which is most easily calculated by replacing the sum over states of discrete mass 

o 47r 2 n 2 



ra- 



ti? 



(3.4) 



by an integral over a continuous mass distribution 5 [9] . The result for time-like propagators 
reads 



D(s) 



„n/2-l 



2M? +2 G N 

and for space-like propagators it is given by 

Uln/2-l 



[vr + 2*J(M a /V5)] 



(3.5) 



with 



I(x) 



D E {t) 



n/2-1 



2M™ +2 G N 



-2i)I E {M a /y/\t\) 



(3.6) 



„2k 



log (x 2 - 1) 



n 



even 



k=l 

(n-l)/2 

£ 

k=l 



(3.7) 



— - — x 2k 1 + - log 
2k - 1 2 & \ x 



x + l 
1 



n = odd . 



5 Note that in principle there is an additional dependence on m| due to the mass-dependent terms in 



Eq. (2.27). However, these mass-dependent terms vanish completely if the particles interacting with the 
graviton are all on-shell. Even if this is not the case, the error resulting from neglecting the additional mj 
dependence in Eq. (3.3) is rather small - our numerical tests showed deviations mostly of the order of the 
numerical precision (~ 1CP 14 ) and in some extreme cases up to 10~ 5 . An effective mass could be chosen in 
such a way that the contribution of the —V term is correct. 
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Ie(x) 



(_l)(n-l)/2 



n/2-1 

E 



(-1)" 



k=l 
'(n-l)/2 

E 

fc=i 



2A: 



+ - log (x 2 + 1) 



(~l) fc 
2fc- 1 



x 2 ^ 1 + tan~ 1 (x) 



n = even 



n = odd . 



(3.8) 



In Eqs. (3.5) and (3.6) above, we have used Eq. ( |3.1| ). In both cases, the leading order 
terms for M 2 3> s are: 



^(1*1 



2MjG N 8 1 ^ 



2M S 4 GW n - 2 



n = 2 



n > 2. 



(3.9) 



By default, we have implemented in the program the complete expressions (3.5 3~8|). 
However, for the sake of comparison with previous work on the subject, we have also im- 
plemented the leading order approximation ( |3.9| ) , as well as the two other scale conventions 
mentioned earlier. In the GRW convention [8], the sum (|3.3| ) amounts to 



D 



while in the Hewett convention [13], it reads 



D 



-iX 



n(MJ H ) 4 G N : 



(3.10) 



(3.11) 



where the unspecified sign A of the amplitude reflects the arbitrariness in the regularization 
procedure employed in obtaining (j3~5|-3i 



3.1.2 External KK particles 

For real graviton production processes, we must sum over all possible final state gravitons. 
A single mass state is labeled by the n-dimensional vector n = (ni,n2, ...,n n ), where the 
rii are integer values. The corresponding mass is given by Eq. ( |3.4|) . The total number 
of accessible graviton mass states in the KK tower is 



TO 



n M™ +2 G N ' 



(3.12) 



with m the maximal kinematically allowed graviton mass. The scale Ms is defined by 
Eq. (|3.1| ). Since this number is enormous, it is impractical to introduce a separate particle 
label in the program for each graviton of definite mass. Instead, it is easier to treat the 
graviton tower as a single particle of variable mass, and then add an extra summation over 
the graviton mass [28,33]. In this respect, real graviton production is unlike any of the 
conventional processes implemented in the general purpose event generators, where the 
particle masses are held fixed. 
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The summation over the graviton mass can be easily carried out by the Monte Carlo 
phase-space integrator used in [59]. Any given scattering event produces a single graviton 
whose mass m^j is picked at random from a uniform distribution of n. Thus, to get the 
total cross section, the weight of every phase space point has to be re-weighted by the 
factor ( gig ). 

3.2 Vertex Feynman Rules 

We have implemented all three- and some of the four-point vertices from [9] for both the 
spin-2 gravitons and the spin-0 dilaton. We also included a triple-Higgs-KK vertex, which 
together with a four-Higgs-KK vertex, was not previously listed in the literature, to the 
best of our knowledge. The Feynman rules for these two additional vertices are shown in 
Figs. and g. 



• 3 m h 
6ia;/c— 



Figure 1: The triplc-Higgs-KK vertex. 




Spin-2 KK: ^ 



Spin-0 KK: -6icoK^J 



Figure 2: The four-Higgs-KK vertex. 



Appendix [A] contains a list of all implemented vertex Feynman rules, translated into 
the functions defined as building blocks in Section 

Note that we define the coupling constants k and u somewhat differently from [9] . This 
should not become an issue, since k and uj are auxiliary variables, and the ADD model is 
defined entirely through M$ and n. In the revised version of [9], the spin-2 propagator and 
the normalization of the polarization tensors have been changed by a factor of 2, which 
has been absorbed into k. We have also rescaled uj by an extra factor of \Jn — 1, which 
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Figure 3: The AMEGIC++ result for e + e~ — > for centre-of-mass energy Ecu = 500 GeV, 

Mg H — 1500 GeV, and A = +1 (crosses) or A = — 1 (squares). The solid histogram is the pure 
Standard Model result. 



typically arises from the multiplicity of the spin-0 KK particles 



K=y/32*G N , W = ^J1_^. (3.13) 

With these definitions the vertex Feynman rules can be simply read off from either the 
published or revised version of [9] and we do not repeat them here. 



4. Tests 



As a test of our implementation we performed a number of checks, comparing AMEGIC++ 
results with independent analytical calculations in the literature: 

• In order to test 2-fermion production processes mediated by virtual graviton ex- 
change, we calculated the angular distributions for e + e~ — > e + e~ — > bb and 
e + e~ — ► cc. We found each one in good agreement with the original work of [13]. 
Fig. |3| shows the AMEGIC++ result for e + e~ — > for a centre-of-mass energy 
Ecm = 500 GeV. We ran AMEGIC++ in the effective operator mode corresponding to 
the Hewett convention and chose Mg = 1500 GeV. Fig. || displays our results for 
both A = +1 (crosses) and A = — 1 (squares), as well as the pure Standard Model 
result (solid histogram). We see perfect agreement with Fig. 1 of the archived ver- 
sion of [13]. We also ran AMEGIC++ in its default mode, choosing the appropriate 
Ms according to Eq. ( p.ll| ) and recovered the curve corresponding to A = +1, thus 
confirming the result of [3]. 
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I I I I I I I I I I I I I I I 1 

500 1000 1500 2000 



E CM (GeV) 

Figure 4: The total cross-section for e + e~ — ► jG as a function of the centre-of-mass energy Eqm, 
for Ms — 2500 GeV, n = 2 and photon energy and angular acceptance cuts of E 1 > 10 GeV and 
| cos# 7 | < 0.9, respectively. Compare with Fig. 3 of [10]. 

• The unpolarized angular distribution for e + e~ — > W + W~ was calculated and found 
to be consistent with the results of [17]. 

• All of our results on 77 — ► hh were found to be in perfect agreement with those of 
Ref. [19]. 

• The total cross sections for 77 — > 77 and e + e~ — > 77 were calculated and compared 
with the results of [16]. We find agreement, except for the case n = 2 in 77 —* 77. We 
believe that the reason for the deviation is that [16] used the approximation \t\ ~ s 
for i-channel diagrams, as in Eq. (p.Sj), while the more general i-channel result is 
given by Eq. (|3.6|). 

• We tested real graviton production through e + e~ — > jG. Our results for the total 
cross section as a function of the centre-of-mass energy Eqm ar e shown in Fig. [| 
for Ms = 2500 GeV, n = 2, and for photon energy and angular acceptance cuts of 
£ 7 > 10 GeV and | cos# 7 | < 0.9, correspondingly. Fig. |4] agrees very well with the 
analogous Fig. 3 of [10] 6 . 

• Finally, we calculated the decay widths of massive KK particles. For spin-2 KK states 
we agree with the results of [9]. For spin-0 particles we agree with the results in the 

6 Note that the calculation in [10] is based on Feynman rules from an earlier version of [9], leading to an 
intermediate result that is too small by a factor of 1/2 [7,66]. However, along with the correction of that 
factor of 1/2, the revised version of [9] also redefined the scale Ms in such a way that the physical cross 
sections for spin-2 KK gravitons remained unchanged (after integration over the KK tower). Thus, Fig. 3 
of [10] can be safely used for comparison purposes here. 
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original printed version of [9], while in its revised version the widths are too large by 
a factor of 2. 

This implementation of the ADD model will be available in the upcoming version 2.0 of 
AMEGIC++. 

5. Conclusions 

The idea that quantum gravity can be relevant already at the TeV scale is extremely at- 
tractive - to theorists and experimentalists alike. Since gravity is universal, there exist 
many potential graviton signatures, but so far only the simplest ones have been explored. 
The calculations become prohibitively complicated and the need for automation soon be- 
comes evident. Thus the main purpose of this paper was to develop the helicity formalism 
for spin-2 particles, which could then be used in an automatic event generator, capable of 
calculating arbitrary processes with gravitons. 

Having completed this first task, we proceeded to implement the ADD model of large 
extra dimensions in AMEGIC++. The large extra dimensions paradigm will soon be tested in 
collider experiments at the Tevatron, LHC and NLC. We believe that having a simulation 
tool like AMEGIC++ at our disposal will greatly facilitate the study of new physics related 
to quantum gravity or extra dimensions. With the general formalism at hand, it is now 
possible to generalize our implementation to models with warped extra dimensions as well. 
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A. Appendix: Vertex Functions 

The following shorthands are used: 

Y = Y(ko, A , feo, A' , c L , c R ) , 

-^oi = X{ko, A , ki, k' , A' , c L , c R ) , 

Xq k = X(k , A , e K , k , A' , cl,c r ) , 

Z 01 = Z(k , \ ,k' ,\' , ki,X 1 ,k' 1 ,X' 1 ,c° L , c R ,c\,c R ) . 

All spin-2 propagators are "cut" (see |2.4| , |2.5| ), while the spin-1 propagators are attached to 
a fermion line. For the case when they are attached to something else in a diagram (and for 
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massive vector bosons), they must be "cut" as well. This is done by simple replacements 
(e.g. to cut the boson line with momentum po): 

Xoi - eX , 
Zqi — ► Xix , 

where A labels the polarization modes of the cut line. 




Figure 5: The fcrmion-fermion-KK vertex. 



FFT = X ox (k -e K - k' -e K ) + X 0ti (k -e x - fc -e A ) , 

-2e x -e K (X(k ,k ,k' ) - X (k ,k' ,k' ) - 2m*Y ) (A.l) 




(A.2) 



Figure 6: The two-gauge-boson-KK vertex. 



VVT = {m\ +p -pi) (x 0X X lK + X 1X X 0K - e A -e K Z i) 

+ (Xoi + X 00 ) (Xio + In ) e A • + Z 01 (p • e K Pi • e A + Pi • e K Po • e A ) 

- (Xoi + X 00 ) po ■ e x + X n po ■ e K ) 

- (X 10 + In) (Xo K pi-e x + XoxW^) , (A.3) 
VVGS = m\Z i + XooX l2 + X 02 Xu . (A.4) 
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/ 
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s 

/ 

/ 

✓ 

/ 

Figure 7: The scalar-scalar-KK vertex. 



55T= ( m | + p . m ) e A . e --p . e A p 1 . e «-p 1 . e A p . e -, (A.5) 
= -2m| + po -Pi • (A.6) 



^ a(pi) Spin-2 KK: {gT^K FFVT 




e "' e » Spin-OKK: —i-gT^uJK FFVGS 



Figure 8: The fermion-fermion-gauge-boson-KK vertex. 



FFVT = X 0K X 1X + X lK X 0X - 2 e A • e K Z 01 , (A.7) 
FFVGS = Z i . (A.8) 
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VVVT = A 012 + A 120 + A 2W 



(A.9) 



A)12 — ( -^OA^lre + -^IA-^Ok ~ e X -€ K Zq\ ) {X 2 § — X 



-21, 



01 



X 2X ( Po -e K - p v e K ) + X 2K ( Po -e x - Pl -e x ) 



Spin-2 KK: -i§«^ e A -e K 




e ^' e ^ Spin-OKK: -Qiujn^ 



Figure 10: The triple-scalar-KK vertex. 
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